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ABSTRACT 

The  ISO  LWS  far-infrared  spectrum  of  the  ultraluminous  galaxy  Mrk  231  shows  OH  and  H20  lines  in  absorption 
from  energy  levels  up  to  300  K  above  the  ground  state,  and  emission  in  the  [O  i]  63  jtm  and  [C  n]  158  jtm  lines.  Our 
analysis  shows  that  OH  and  H20  are  radiatively  pumped  by  the  far-infrared  continuum  emission  of  the  galaxy.  The  ab¬ 
sorptions  in  the  high-excitation  lines  require  high  far-infrared  radiation  densities,  allowing  us  to  constrain  the  properties 
of  the  underlying  continuum  source.  The  bulk  of  the  far-infrared  continuum  arises  from  a  warm  (Tdust  =  70-100  K), 
optically  thick  (rioo/jm  =  1  -2)  medium  of  effective  diameter  200-400  pc.  In  our  best-fit  model  of  total  luminosity  Fir, 
the  observed  OH  and  H20  high-lying  lines  arise  from  a  luminous  {L/L  ir  ~  0.56)  region  with  radius  ~  1 00  pc.  The  high 
surface  brightness  of  this  component  suggests  that  its  infrared  emission  is  dominated  by  the  AGN.  The  derived  column 
densities  N( OH)  >  1 0 1 7  cm~2  and  N(  H20)  >6x10' 6  cm~2  may  indicate  X-ray  dominated  region  (XDR)  chemistry, 
although  significant  starburst  chemistry  cannot  be  ruled  out.  The  lower-lying  OH,  [C  n]  158  jtm,  and  [O  i]  63  /im  lines 
arise  from  a  more  extended  (~350  pc)  starburst  region.  We  show  that  the  [C  n]  deficit  in  Mrk  23 1  is  compatible  with  a 
high  average  abundance  of  C+  because  of  an  extreme  overall  luminosity  to  gas  mass  ratio.  Therefore,  a  [C  n]  deficit  may 
indicate  a  significant  contribution  to  the  luminosity  by  an  AGN,  and/or  by  extremely  efficient  star  formation. 

Subject  headings:  galaxies:  abundances  —  galaxies:  individual  (Mrk  231)  —  galaxies:  ISM  — 
galaxies:  starburst  —  infrared:  galaxies  —  radiative  transfer 


1.  INTRODUCTION 

The  peculiar  ultraluminous  infrared  galaxy  (ULIRG;  Fir  > 
1012  Lr.j)  Markarian  231  (Mrk  231,  12540+5708)  is  the  most 
luminous  infrared  galaxy  in  the  local  universe,  with  a  8- 1 000  jtm 
luminosity  of  3.2  x  10 12  Le  (Sanders  et  al.  2003),  and  may  be  a 
representative  example  of  the  link  between  active  galactic  nuclei 
(AGNs)  and  nuclear  starbursts  (Scoville  2004).  A  QSO-like  nu¬ 
cleus  is  evident  from  many  observations:  optically  it  is  classified 
as  a  Type  1  Seyfert  (Boksenberg  et  al.  1977;  Cutri  et  al.  1984; 
Baan  et  al.  1998),  it  exhibits  UV  through  IR  polarization  and  broad 
absorption  lines  (Smith  et  al.  1995),  it  has  compact  X-ray  emission 
(e.g.,  Gallagher  et  al.  2002)  and  extremely  compact  mid-infrared 
emission  (Soifer  et  al.  2000),  and  in  the  radio  it  is  variable  and 
possesses  a  parsec-scale  jet  (Ulvestad  et  al.  1999;  Taylor  et  al. 
1999).  Nevertheless,  there  is  also  evidence  of  a  compact  starburst 

1  Based  on  observations  with  the  Infrared  Space  Observatory,  an  ESA  project 
with  instruments  funded  by  ESA  Member  States  (especially  the  principal  investi¬ 
gator  countries:  France,  Germany,  Netherlands,  and  the  United  Kingdom)  and  with 
the  participation  of  ISAS  and  NASA. 


in  these  results  as  well  as  in  Very  Large  Array  (VLA)  observa¬ 
tions  of  H  i  21  cm  absorption  (Carilli  et  al.  1998),  near-infrared 
observations  ( Tacconi  et  al.  2002),  and  millimeter  CO  interfer¬ 
ometry  (Bryant  &  Scoville  1996;  Downes  &  Solomon  1998, 
hereafter  DS98).  Estimates  for  the  starburst  luminosity  range 
from  1/3  to  2/3  of  the  bolometric  luminosity  (Davies  et  al.  2004; 
DS98). 

Molecular  observations  have  provided  important  clues  about 
the  concentration  and  kinematics  of  the  gas  in  Mrk  231.  DS98 
showed  the  presence  of  an  inner  nuclear  disk  of  radius  ^460  pc 
in  CO(2- 1 ),  and  a  more  extended  disk  with  lower  brightness.  Most 
of  the  molecular  gas  has  been  found  to  be  dense  (~  l  04  cm~3) 
and  warm  (^70  K)  from  recent  observations  of  CO  and  HCN  sub- 
millimeter  lines  (Papadopoulos  et  al.  2007,  hereafter  PIW07). 
Lahuis  et  al.  (2007)  have  inferred  embedded  starburst  chemistry 
in  Mrk  23 1  and  other  ULIRGs  based  on  mid-lR  Spitzer  observa¬ 
tions  of  rovibrational  bands  of  warm/hot  HCN  and  C2H2,  while 
Gracia-Carpio  et  al.  (2006)  and  Aalto  et  al.  (2007)  have  inferred 
XDR  chemistry  and/or  radiative  pumping  based  on  anomalous 
intensity  ratios  of  millimeter  lines  of  HCN,  HNC,  and  HCO+. 
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The  bulk  of  the  luminosity  in  ULIRGs  is  emitted  at  far- infrared 
(FIR)  wavelengths,  where  a  number  of  molecular  tracers  are  de¬ 
tected,  mostly  in  absorption.  Prominent  lines  of  OH  and  H20  were 
detected  using  Infrared  Space  Observatory  (ISO)  LWS  in  the  FIR 
spectrum  of  Arp  220,  along  with  absorption  features  by  radicals 
such  as  NH  and  CH,  revealing  a  chemistry  that  may  be  indicative 
of  photodissociation  regions  ( PDRs)  with  plausible  contribution 
by  shocks  and  hot  cores  (Gonzalez- Alfonso  et  al.  2004,  hereafter 
Paper  I).  However,  those  species  are  also  expected  to  be  enhanced 
in  XDRs  (Meijerink  &  Spaans  2005),  so  that  the  dominant  chem¬ 
istry  in  the  nuclear  regions  of  ULIRGs  remains  uncertain.  In 
Paper  I,  the  ISO/ LWS  FIR  spectrum  of  Arp  220  was  analyzed 
by  means  of  radiative  transfer  calculations,  which  included  a 
nonlocal  treatment  of  the  molecular  excitation  by  absorption  of 
FIRphotons.  Paper  I  showed  that  the  population  of  high-excitation 
OH  and  H20  rotational  levels,  in  evidence  from  absorption  in 
high-lying  lines,  is  pumped  through  absorption  of  FIR  continuum 
photons,  a  process  that  requires  high  FIR  radiation  densities.  The 
detection  of  these  lines  thus  not  only  reveals  the  chemical  and 
excitation  conditions  in  the  absorbing  regions,  it  also  sheds  light 
on  the  size  and  characteristics  of  the  underlying  continuum  FIR 
source  in  spite  of  the  low  angular  resolution  currently  available  at 
these  wavelengths. 

In  this  paper  we  extend  our  approach  of  Paper  I  to  the  ISO/  LW S 
FIR  spectrum  of  Mrk  231,  and  show  that  this  galaxy  spectrum 
presents  striking  similarities  to  that  of  Arp  220.  Specifically,  strong 
absorption  in  the  high-excitation  OH  and  H20  lines  is  also  seen 
in  Mrk  23 1 .  Rotationally  excited  OH  in  Mrk  23 1  has  been  pre¬ 
viously  detected  via  the  21 1 1/2  A-doublet  transitions  (Henkel  et  al. 
1987).  Very  long  baseline  interferometry  ( VLB1)  observations  of 
the  megamaser  OH  emission  at  1 8  cm  wavelength  trace  an  inner 
torus  or  disk  of  size  ~  100  pc  around  the  AGN  (Klockner  et  al. 
2003),  and  MERLIN  observations  were  able  to  map  essentially 
the  whole  single-dish  megamaser  OH  emission  with  angular  reso¬ 
lution  of  «0.3"  (Richards  et  al.  2005).  We  analyze  here  both  the 
FIR  continuum  emission  and  the  high-excitation  OH  and  H20 
lines,  as  well  as  the  [C  n]  158  pm  and  [O  i]  63  pm  emission  lines. 
In  §  2  we  present  the  ISO  spectroscopic  observations  of  Mrk  23 1 . 
In  §  3  we  first  analyze  simple  models  for  the  FIR  continuum  emis¬ 
sion  from  Mrk  231,  and  then  examine  how  well  those  models 
reproduce  the  observed  FIR  emission  and  absorption  lines.  Sec¬ 
tion  4  summarizes  our  results.  We  adopt  a  distance  to  Mrk  23 1  of 
170  Mpc  (H0  =  75  km  s_1  Mpc~]  andr  «  0.042). 

2.  OBSERVATIONS  AND  RESULTS 

The  frill  43-197  pm  spectrum  of  Mrk  231  (first  shown  and 
discussed  by  Harvey  et  al.  1999),  was  obtained  with  the  LWS 
spectrometer  (Clegg  etal.  1996)  on  board  ISO  ( Kessler  etal.  1996). 
In  Figure  1 ,  it  is  compared  with  that  of  Arp  220  (Paper  I)  rescaled 
to  the  same  distance  (170  Mpc).  The  grating  spectral  resolution  is 
~0.3  pm  in  the  43-93  pm  interval  (detectors  SW1-SW5),  and 
~0.6  pm  in  the  80-197  pm  interval  (detectors  LW1-LW5),  cor¬ 
responding  to  Av  ^  103  km  s  1 .  The  lines  are  thus  unresolved  in 
velocity  space.  The  «80"  beam  size  ensures  that  all  the  FIR  con¬ 
tinuum  and  line  emission/absorption  from  Mrk  23 1  (CO  size  ~4"; 
DS98)  lies  within  the  ISO  LWS  aperture. 

The  data  (TDT  numbers  5100540, 18001306,  and  60300241) 
were  taken  from  the  highly  processed  data  product  (HPDP)  data 
set  (called  uniformly  processed  LWS01  data),  and  reduced  using 
version  10.1  of  the  Off  Line  Processing  (OLP)  Pipeline  system 
(Swinyard  et  al.  1996).  We  performed  subsequent  data  process¬ 
ing,  including  co-addition,  scaling,  and  baseline  removing,  using 
the  ISO  Spectral  Analysis  Package  (IS  AP;  Sturm  etal.  1998)and 
our  own  routines.  In  order  to  obtain  a  smooth  spectrum  through- 
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Fig.  1 . — Comparison  between  the  FIR  emission  from  Mrk  23 1  {black  line ) 
and  Arp  220  {gray  line).  The  spectrum  of  Arp  220  has  been  rescaled  to  the  dis¬ 
tance  of  Mrk  23 1  (170  Mpc).  The  position  of  the  OH  lines  discussed  in  the  text 
are  indicated  and  labeled  with  their  wavelengths. 


out  the  whole  LWS  range,  the  flux  densities  given  by  each  de¬ 
tector  were  corrected  by  multiplicative  scale  factors.  Corrections 
were  less  than  25%  except  for  detectors  LW2  and  LW3  ( 1 00— 
145  pm),  for  which  the  corrections  were  30%.  We  thus  attribute 
an  uncertainty  of  30%  to  the  overall  continuum  level,  as  well  as 
for  the  line  fluxes. 

Figure  1  shows  that  the  FIR  spectra  of  Mrk  23 1  and  Arp  220 
are  similar  in  key  aspects  (see  also  Fischer  et  al.  1 999),  in  particular 
the  prominent  molecular  absorptions  mostly  due  to  OH  doublets 
(that  will  be  referred  to  hereafter  as  lines)  and  the  lack  of  strong 
fine-structure  line  emission  typically  seen  in  less  luminous  galaxies. 
A  closer  inspection  of  the  pattern  of  line  emission/absorption 
in  both  sources  is  shown  in  Figure  2,  where  the  continuum- 
normalized  spectra  are  compared.  Of  particular  interest  are 
the  clear  detections  in  both  sources  of  the  high-excitation  OH 
n3/27/2-5/2  84  pm  and  n3/?9/2-7/2  65  pm  lines,  with  lower 
level  energies  of  120  and  290  K,  respectively  (see  §  3).  The 
330  — >  22i  and  331  — >  2?o  H20  66-67  pm  lines,  both  with  lower 
levels  at  195  K,  are  also  detected  in  Mrk  23 1,  as  well  as  the  ten¬ 
tatively  identified  22o  — >  In  line  at  101  pm.  It  is  likely  that  the 
increased  noise  level  at  A  <:  160  pm  is  responsible  for  the  non¬ 
detection  of  the  high-excitation  Hi/23/2 — 1/2  OH  line  in  Mrk 
231,  which  is  seen  in  strong  emission  in  Arp  220.  While  the  high- 
excitation  OH  and  H20  lines  at  65-67  pm  are  of  similar  strength 
in  Mrk  23 1  and  Arp  220,  the  H20  lines  at  longer  wavelengths 
are  undoubtedly  weaker  in  Mrk  231,  as  seen  for  the  322  — >  2n, 
220  — ^  In,  and  221  — >  1 10  H20  lines  at  90,  102,  and  108  pm, 
respectively.  The  weakness  of  the  latter  lines  in  Mrk  23 1  suggests 
that  the  region  where  the  high-lying  H20  lines  are  formed  is  rela¬ 
tively  weak  in  the  FIR  continuum  at  X  =  90-108  pm.  The  Mrk  23 1 
spectrum  thus  suggests  that  a  warm  component  with  relatively 
weak  contribution  to  the  FIR  continuum  at  X  S;  80  pm  is  respon¬ 
sible  for  the  observed  high-excitation  absorptions  (§  3.3).  Table  1 
lists  the  line  fluxes,  continuum  flux  densities  at  the  correspond¬ 
ing  wavelengths,  and  equivalent  widths  for  the  lines  detected  in 
Mrk  231. 
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Fig.  2. — Comparison  between  the  continuum-normalized  spectra  of  Mrk  23 1  ( upper  spectra,  black  lines)  and  Arp  220  {lower  spectra,  gray  lines).  The  positions  of 
lines  discussed  in  the  text  are  indicated  and  labeled. 


In  the  case  of  Arp  220,  we  used  high  spatial  resolution  con¬ 
tinuum  measurements  available  in  the  literature  to  infer  that 
Arp  220  is  optically  thick  even  in  the  submillimeter  continuum 
(Paper  I;  see  also  Downes  &  Eckart  2007).  The  steeper  decrease 
of  the  flux  density  with  increasing  wavelength  in  Mrk  23 1 ,  how¬ 
ever,  suggests  that  it  has  lower  FIR  continuum  opacities  (Fig.  1). 
This  expectation  is  further  reinforced  by  the  detection  in  Mrk  23 1 
of  the  [N  n]  122  /./m  line,  a  feature  not  seen  in  Arp  220  (Fig.  2). 
Other  notable  differences  between  both  sources  are  that  the  [O  i] 
63  pm  line  is  observed  in  emission  in  Mrk  23 1  but  in  absorption  in 
Arp  220,  and  that  the  ground-state  119,  53,  and  79  /im  OH  lines 
are  significantly  weaker  in  Mrk  23 1  (Fig.  2).  In  modeling  Arp  220, 
we  were  forced  to  invoke  an  absorbing  “halo”  to  account  for 
these  lines;  in  Mrk  231,  no  such  halo  is  required  (§  3). 

In  the  spectrum  of  Mrk  231,  the  main  119.3  pm  OH  line  ap¬ 
pears  to  be  slightly  blueshifted  relative  to  the  expected  position, 
an  effect  we  attribute  to  the  proximity  of  the  line  to  the  edge  of 
the  LW3  detector.  There  is  a  nearby  weaker  redshifted  feature 
at  120  /im  which  coincides  with  the  expected  position  of  the 
ground  n  3/2  5/2 -3/2  I8OH  line,  and  appears  as  a  marginal  fea¬ 
ture  in  both  the  “up”  and  “down”  grating  scans.  However,  the 
limited  signal-to-noise  ratio  [(1.0  ±  0.4)  x  10~20  W  cm-2],  the 
narrow  appearance  of  the  feature  («0.42  /im),  and  the  fact  that  it 
is  not  blueshifted  as  the  main  line,  make  that  assignment  only 
tentative.  In  Arp  220,  the  main  OH  line  is  not  shifted  because  it 


does  not  fall  so  close  to  the  edge  of  the  detector,  as  a  conse¬ 
quence  of  the  lower  redshift  of  the  source.  In  Arp  220,  a  redshifted 
shoulder  appears  at  120  /im,  suggesting  the  possibility  that  18OH 
may  be  responsible  for  it  (Paper  1).  We  cannot  however  be  certain 
that  18OH  is  detected  in  any  of  these  sources,  but  given  the  high 
16OH  column  densities  we  derive  in  some  of  our  models  below 
(§  3.3)  and  the  fact  that  values  of  the  16OH/18OH  ratios  below 
the  canonical  value  of  500  may  be  expected  in  regions  where 
the  interstellar  matter  (ISM)  is  highly  processed  by  starbursts 
(Paper  I),  our  tentative  identification  should  be  followed  up  with 
future  Herschel  Space  Observatoiy  observations  with  higher  spec¬ 
tral  resolution  and  sensitivity.  Finally,  the  spectrum  of  Mrk  231 
shows  a  broad  feature  at  the  position  of  the  ni/2-n3/23/2-3/2 
OH  line  (53  /im).  We  note  that  the  blueshifted  side  of  this  ab¬ 
sorption  is  coincident  with  the  OH  n 3/2 11/2-9/2  line,  with  a 
lower  level  energy  of  5 1 1  K;  however,  the  proximity  of  this  spec¬ 
tral  feature  to  the  edge  of  the  SW2  detector  precludes  any  defin¬ 
itive  assignment. 

The  FIR  detections  of  both  NH  and  NH3  in  Arp  220  were 
reported  in  Paper  I.  NH3  was  also  detected  via  the  25  GHz  in¬ 
version  transitions  by  Takano  et  al.  (2005),  who  derived  a  NH3 
column  density  6  times  higher  than  our  value.  The  difference  likely 
arises  because  of  the  high  FIR  continuum  opacities  in  Arp  220, 
which  cause  the  observed  FIR  absorptions  to  trace  only  a  fraction 
of  the  total  gas  column.  Since  there  are  no  such  extinction  effects 
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TABLE  1 

Line  Fluxes,  Continuum  Flux  Densities,  and  Equivalent  Widths  for  the  Lines  Detected 
in  the  ISO  LWS  Spectrum  of  Mrk  231 


2resta  Line  Fluxb  Continuum0  WA 

Species  Transition  (/.on)  (10'2°  W  cm-2)  (10'19  W  cm_2fim_1)  (10_2/tm) 

OH .  ni/2-n3/2  3/2-3/2  53.3  -7.0  ±  2.0  33.5  2.1  ±  0.6 

[0  i] .  7’  iP2  63.2  2.6  ±  0.8  25.0  -1.0  ±  0.3 

OH .  II,.,  D  2  7  2  65.2  -4.1  ±  0.6  23.8  1.7  ±  0.3 

H20  .  330-221  66.4  -3.3  ±  0.9  23.1  1.4  ±  0.4 

H20  .  33i-220  67.1  -2.2  ±  0.9  22.7  1.0  ±  0.4 

OH .  I  Is  2  7/2— 5/2  84.5  -8.6  ±  0.7  13.8  6.2  ±  0.5 

H20 .  220—  1 1 1  101  -1.3  ±  0.3  8.70  1.5  ±  0.4 

OH .  n3/2  5/2— 3/2  119.3  -3.6  ±  0.3  5.22  6.9  ±  0.6 

[Nn] .  3T>2-3Pi  121.8  1.5  ±  0.2  4.72  -3.2  ±  0.4 

[C  II] .  2Pin-2Pm  157.7  3.7  ±  0.1  1.73  -21.4  ±  0.6 


a  For  OH  doublets,  an  average  for  the  two  components  is  given. 

b  Errors  do  not  include  calibration  uncertainties  in  the  continuum  level.  Negative  (positive)  values  indicate  that  the  line  is 
detected  in  absorption  (emission). 

c  Uncertainties  in  the  continuum  level  are  less  than  30%. 

d  Equivalent  widths  are  positive  (negative)  for  lines  observed  in  absorption  (emission). 


at  25  GHz,  the  NH3  inversion  transitions  are  expected  to  trace 
higher  NH3  column  densities.  Figure  2  shows  that,  by  contrast, 
the  NH3  lines  are  not  detected  in  Mrk  23 1 ,  although  the  relatively 
high  noise  at  125  /im  does  not  rule  out  future  detection  of  NH3 
with  Herschel  at  a  level  similar  to  that  of  Arp  220. 

There  are  two  marginally  detected  (2.5  a  level )  spectral  features 
seen  at  153.0  and  152.3  /im  in  the  Mrk  231  spectrum  (Fig.  3).  Al¬ 
though  close  to  the  expected  position  of  the  main  NH  feature  at 
153.22  /im,  the  153.0  /on  feature  appears  significantly  shifted  by 
0.25  /im  from  it,  and  better  coincides  with  the  position  of  the 
OH+23-l2  line.  Also,  the  152.3  /tm  feature  lies  at  0.1  /im  from 
the  expected  position  of  the  OH+22-l  i  line.  In  Paper  1,  we  also 
suggested  that  OH+  could  contribute  to  the  spectrum  of  Arp  220 
for  two  reasons:  ( 1 )  our  models  were  unable  to  reproduce,  using 
NH  and  NH3,  the  observed  strong  absorption  at  102  /im,  which 
coincides  with  the  expected  position  of  the  OH+34-23  line; 
(2)  there  was  an  absorption  feature  at  76.4  /im  that,  if  real,  could 
be  attributed  to  the  OH+44-33  transition.  Since  OH+  has  never 
been  detected  in  the  Galactic  interstellar  medium  or  that  of  any 
galaxy,  here  we  only  highlight  the  intriguing  possibility  of  its  de¬ 
tection  in  two  ULIRGs.  Sensitive,  higher-resolution  Herschel  ob¬ 
servations  are  needed  to  resolve  this  tantalizing  speculation. 
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Fig.  3. — Comparison  between  the  spectra  of  Mrk  231  {upper  spectrum,  histo¬ 
gram)  and  Arp  220  {lower  spectrum,  solid  line )  around  152  pm.  The  spectrum  of 
Arp  220  has  been  rescaled  to  the  distance  of  Mrk  23 1  (1 70  Mpc).  The  position  of  the 
NH,  CH,  and  OH+  lines  are  indicated  and  labeled. 


The  luminosity  of  the  [C  n]  Pyi-2P  M2  fine-structure  line  at 
158  /im  is  2.5  times  stronger  in  Mrk  231  than  in  Arp  220,  but 
given  the  higher  FIR  luminosity  of  this  source  (Fig.  1 ),  the  [C  n] 
to  FIR  luminosity  ratios  are  rather  similar,  with  values  of  2.5  x 
10~4  and  2.1  x  10~4  for  Mrk  231  and  Arp  220,  respectively 
(Luhman  et  al.  2003).  These  are  among  the  lowest  values  found 
in  galaxies,  illustrating  the  so-called  [C  n]  deficit  found  in  ULIRGs. 
The  [C  n]  line  emission  from  Mrk  231  is  analyzed  in  §  3.4. 

3.  ANALYSIS 

3.1.  Models  for  the  FIR  Continuum 

Figure  4  illustrates  several  ways  that  the  FIR  to  millimeter 
continuum  can  be  fit  and  interpreted.  We  first  modeled  (model  A 
in  Fig.  4a)  the  FIR  source  in  Mrk  23 1  as  an  ensemble  of  identical 
dust  clouds,  each  of  which  is  heated  by  its  own  single  central 
luminosity  source.  The  representative  cloud  is  assumed  to  be  spher¬ 
ical,  with  radius  Rc,  and  is  divided  into  concentric  shells  whose 
dust  temperatures  are  computed  from  the  balance  of  heating  and 
cooling  (Gonzalez- Alfonso  &  Cemicharo  1999).  We  used  a  mix¬ 
ture  of  silicate  and  amorphous  carbon  grains  with  optical  con¬ 
stants  from  Preibisch  et  al.  (1993)  and  Draine  (1985).  The  stellar 
continuum  was  taken  from  Leitherer  et  al.  (1999),  but  results 
depend  only  weakly  on  this  choice  because  the  intrinsic  contin¬ 
uum  is  absorbed  by  the  dust  and  re-emitted  at  infrared  wave¬ 
lengths.  Once  the  equilibrium  temperatures  are  obtained  for  each 
shell,  the  resulting  continuum  emission  from  the  cloud  is  com¬ 
puted  and  multiplied  by  Nc,  the  number  of  clouds  in  the  source 
required  to  match  the  absolute  flux  densities.  This  scaled  spec¬ 
trum  is  shown  in  Figure  4 a.  The  other  three  models  (B,  C,  and  D, 
shown  in  Figs.  4b- 4d)  use  graybodies  with  uniform  dust  tem¬ 
peratures  Td  to  characterize  the  continuum  emission  (e.g.,  Roche 
&  Chandler  1993;  Annus  et  al.  2007). 

Assuming  that  the  individual  clouds  do  not  overlap  along  the 
line  of  sight,  our  results  do  not  depend  particularly  on  the  radius 
or  luminosity  adopted  for  the  model  individual  cloud  because 
identical  results  are  obtained  if  Rc  is  multiplied  by  a  factor  of  a, 
the  luminosity  by  a2,  Nc  by  a~2,  and  the  continuum  opacity  is 
kept  constant  (see  Paper  1).  The  models  are  thus  characterized 
by  the  luminosity  of  the  whole  ensemble,  the  radial  opacity  of 
the  clouds  at  a  given  wavelength  (which  we  adopt  to  be  100  /im: 
t loo/im)?  and  the  equivalent  radius  of  the  source,  defined  as  Re0i  = 
NpRc.  These  parameters  are  listed  in  Table  2. 
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Fig.  4. — Different  models  for  the  FIR  continuum  emission  from  Mrk  231.  The  whole  set  of  parameters  that  characterize  models  A,  B,  C,  and  D  are  listed  in  Table  2.  Flux 
densities  at  800  and  1100  fin i  are  taken  from  Roche  &  Chandler  (1993;  corrected  for  nonthermal  emission),  at  450  /im  from  Rigopoulou  et  al.  (1996),  at  350  fim  from  Yang  & 
Phillips  (2007),  and  at  2  <  25  fim  from  Rieke  (1976). 

In  model  A,  the  individual  clouds  are  optically  thin  so  that  temperature  through  the  whole  region  (disk),  so  that  the  resulting 
some  degree  of  cloud  overlap  would  yield  a  similar  fit  to  the  extent  of  the  FIR  emission  is  Req  =  NlnRc. 

continuum  while  decreasing  the  value  of  Xeq-  For  instance,  if  the  The  observed  continuum  can  be  reproduced  from  model  A’s 

clouds  are  distributed  in  a  spherical  volume,  Req  =  7VC1/3RC,  giving  cloud  ensemble  that  is  optically  thin  in  the  FIR.  Model  A  also 

Req  =  400  pc  for  clouds  with  Rc  =  20  pc.  However,  the  predicted  predicts  that  the  starburst  dominates  the  continuum  for  A  S;  15/im, 

opacity  through  the  modeled  region,  yV(!/3rioo/im,  will  be  much  while  the  torus/disk  around  the  AGN  would  then  dominate  the 

higher  than  that  of  an  individual  cloud,  and  this  physical  situation  mid-infrared  continuum,  in  qualitative  agreement  with  the  mod- 

is  already  described  in  models  B-C,  where  higher  opacities  along  els  by  Farrah  et  al.  (2003).  The  equivalent  radius  of  the  starburst 

the  line  of  sight  and  a  more  compact  region  of  FIR  emission  are  is  slightly  larger  than  the  radius  of  the  outer  disk  observed  by 

assumed.  In  order  to  avoid  this  model  redundancy,  we  choose  DS98.  Because  riooMm  is  low  andReq  is  high,  this  model  predicts 

our  continuum  models  such  that  an  individual  cloud  describes  that  the  FIR  radiation  density >  is  low,  a  prediction  that  is  not  con- 

the  characteristic  continuum  opacity  (rioo/jm  in  Table  2)  and  dust  sistent  with  our  models  of  the  observed  OH  line  strengths  (§  3.2). 

TABLE  2 

Models  for  the  Continuum  Emission 


Rcq  W(H,)a  Td  rrad  M  L 

Model  Component  Tm,m  (pc)  (cm-2)  (K)  (K)  (M0)  (L0)  i/iIR 

A .  ...  6.8  x  10'2  1  8  3  0  4.6  x  1022  Variable15  23  1.0  x  1010  2.1  •  10  2  0.72 

B .  ...  0.55  4  0  0  3.7  x  1023  5  5  4  4  4.0  x  109  1.7  x  1012  0.50 

C .  ...  2.2  2  00  1.5  x  1024  7  4  7  1  4.0  x  109  2.0  x  1012  0.61 

D .  Warm  1.1  106  7.2  x  1023  1  00  8  5  5.5  x  10s  1.8  x  1012  0.56 

D .  Cold  0.9  7  3  5  0  6.5  x  1023  4  7  4  2  5.3  x  109  7.4  x  10  1  0.22 


A  gas-to-dust  mass  ratio  of  100  is  assumed,  together  with  a  mass-absorption  coefficient  for  dust  of  44  cm2  g  1  at  100  /an. 
The  temperature  is  calculated  from  the  balance  between  heating  and  cooling;  see  text  for  details. 
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As  both  Tioo/im  and  Ta  are  increased  in  models  B  and  C,  the 
radiation  density  increases  and,  therefore,  the  equivalent  size  re¬ 
quired  to  reproduce  the  observed  emission  becomes  smaller.  As 
a  consequence,  models  B  and  C  predict  increasing  compactness 
of  the  dust  clouds  responsible  for  the  FIR  emission,  with  Rsq  = 
400  and  200  pc,  respectively.  With  a  single-component  model, 
however,  Req  cannot  be  reduced  more  than  in  model  C  without 
degrading  the  quality  of  the  fit.  However,  a  two-component  model 
as  shown  in  model  D  is  able  to  reproduce  the  FIR  emission,  in¬ 
voking  a  quite  compact  (~100  pc)  and  warm  ( 100  K)  component 
(f^warm),  and  a  colder  and  more  extended  one  that  dominates  at 
X  >  80  /tm  CD Coid). 

A  convenient  way  to  characterize  the  radiation  density  in 
the  modeled  regions  is  to  compute  the  radiation  temperature  at 
100  fim  from 

his 

Fmd (100  /m n)  =  +  ’  ( !) 

where  f l  =  ttNcR2/D2  is  the  solid  angle  subtended  by  the  mod¬ 
eled  source,  Fioo^m  is  the  predicted  flux  density  at  100  /tm,  and 
other  symbols  have  their  usual  meanings.  7^(100  /im)  is  also 
listed  in  Table  2,  together  with  the  gas  mass,  luminosity,  and  frac¬ 
tion  of  the  bolometric  luminosity  for  each  model.  The  calculated 
gas  masses  assume  a  gas-to-dust  mass  ratio  of  100.  In  all  cases, 
they  are  lower  than  the  dynamical  masses  determined  by  DS98 
when  Req  is  identified  with  the  radial  extent  of  the  source  (and 
therefore  compatible  with  the  inferred  rotation  velocities  in  the 
disk).  Our  inferred  masses  are  consistent  in  models  B-D  with  the 
mass  inferred  by  P1W07,  but  are  in  all  cases  higher,  by  at  least  a 
factor  of  2,  than  the  gas  masses  obtained  by  DS98.  This  dis¬ 
crepancy  may  be  explained  in  at  least  four  possible,  different  ways: 
(1)  the  physical  radial  extent  of  the  cloud  ensemble,  which  ac¬ 
counts  for  cloud  filling,  is  given  by  RT  =  f~l/1Req,  where /is 
the  area  filling  factor,  so  that  Req  is  a  lower  limit  of  Rt;  (2)  our 
calculated  masses  depend  on  the  mass-absorption  coefficient  for 
dust,  which  we  have  assumed  to  be  k  1300pm  =  0.33  cm2  g~ 1  based 
on  a  mixture  of  silicate  and  amorphous  carbon  grains  (Preibisch 
et  al.  1 993 ;  Draine  1985),  but  could  be  up  to  a  factor  ~6  higher  if 
the  dust  is  mainly  composed  of  fluffy  aggregates  ( Kruegel  & 
Siebenmorgen  1994);  (3)  the  gas-to-dust  mass  ratio  may  depart 
significantly  from  the  standard  value  of  100;  and  (4)  the  masses 
derived  by  DS98  for  Mrk  23 1  could  be  lower  limits  in  the  light 
of  the  submillimeter  CO  emission  reported  by  PIW07.  A  com¬ 
bination  of  these  factors  may  explain  our  higher  values. 

The  luminosities  in  Table  2  account  for  50%-80%  of  the  ob¬ 
served  8-1000  /im  infrared  luminosity.  Model  A  implicitly  as¬ 
sumes  that  the  calculated  luminosity  has  a  starburst  origin;  the 
luminosity  from  model  B  and  from  the  cold  component  of  model 
D  are  also  attributable  to  the  starburst  in  view  of  the  spatial  extent 
of  the  modeled  source.  Since  model  C  and  the  warm  compo¬ 
nent  of  model  D  are  more  compact,  a  combination  of  AGN  and 
starburst  contributions  is  more  plausible.  The  surface  brightness 
in  model  C  is  4  x  10 12  Le  kpc~2,  a  factor  of  2  higher  than  the 
peak  global  value  found  in  starburst  galaxies  by  Meurer  et  al. 
(1997),  suggesting  an  important  (but  uncertain)  contribution  by 
the  AGN  to  the  observed  FIR  emission  (Soifer  et  al.  2000).  In 
addition,  the  luminosity-to-mass  ratio  of  500  Le/M@  coincides 
with  the  uppermost  limit  proposed  by  Scoville  (2004)  for  a  star- 
burst.  The  very  high  surface  brightness  (1.3  x  10 13  Le  kpc~2) 
and  luminosity-to-mass  ratio  (^3300  Le/M@ )  of  the  warm  com¬ 
ponent  of  model  D  (//vaml),  as  well  as  its  compactness,  persua¬ 
sively  indicate  that  this  component  is  most  probably  dominated  by 


the  AGN.  The  most  plausible  relative  contributions  by  the  AGN 
and  the  starburst  to  Dw arm  are  discussed  in  §  4. 

In  summary,  different  approaches  can  be  used  to  successfully 
fit  the  observed  FIR  continuum  emission,  with  the  properties  of 
the  clouds  that  emit  that  radiation  in  these  approaches  spanning  a 
wide  range  of  possible  physical  scenarios.  But  ISO  LWS  has 
provided  us  with  spectroscopic  information,  and  we  show  next 
how  the  observed  high  excitation  OH  and  H2O  lines  impose 
important  constraints  on  these  continuum  models. 

3.2.  Equivalent  Widths 

We  analyze  the  OH  equivalent  widths  assuming  that  the  OH 
molecules  form  a  screen  in  front  of  the  IR  source.  The  strengths 
of  the  n3/27/2-5/2  and  9/2— 7/2  OH  doublets  at  84  and  65  /im 
enable  us  to  conclude  that  the  excited  OH  covers  a  substantial 
fraction  of  the  FIR  emission  region.  Assuming  that  each  line  of 
the  84  /im  doublet  absorbs  all  the  background  84  /jin  continuum 
over  a  velocity  range  of 250  km  s_1  along  each  line  of  sight,  and 
that  there  is  no  significant  re-emission  in  the  line,  the  covering 
factor  is  ~50%.  This  value  may  be  considered  a  lower  limit  for 
the  following  reasons.  The  submillimeter  CO  line  profiles  shown 
by  PIW07  have  FWHMs  of  200-250  km  s-1,  and  the  lines  are 
expected  to  be  broadened  by  velocity  gradients  and,  in  particular, 
by  the  disk  rotation;  therefore,  the  velocity  range  of  250  km  s_1 
assumed  above  is  probably  an  upper  limit.  DS98  inferred  local 
turbulent  velocities  of  up  to  60  km  s_1  at  inner  radii  (100  pc)  and 
decreasing  as  r~0'3.  If  we  adopt  an  intrinsic  Gaussian  line  profile 
with  the  highest  value  of  the  turbulent  velocity,  A  V  =  60  km  s_1, 
and  saturate  the  84  /im  line  to  the  degree  that  an  effective  width2  of 
250  km  s~*  is  obtained  for  each  component  of  the  doublet,  the 
derived  84  /t m  foreground  opacity  at  line  center  is  ~50,  but  the 
high  column  density  required  for  this  opacity  is  hard  to  reconcile 
with  that  inferred  from  the  other  observed  OH  line  strengths 
(§  3.3).  Finally,  some  significant  re-emission  in  the  84  /tm  OH 
line  is  expected  because  the  H3/29/2-7/2  OH  line  at  65  /im  that 
originates  from  its  upper  level  is  detected  in  absorption.  We 
therefore  conclude  that  the  observed  84  /tm  OH  absorption  is 
widespread,  and  probably  covers  the  bulk  of  the  84  /tm  continuum 
emission  regions.  On  the  other  hand,  the  opacities  in  the  high- 
lying  65  /tm  line  should  only  be  moderate;  for  reference,  if  we 
adopt  for  each  component  an  upper  limit  of  150  km  s_1  on  the 
effective  velocity  interval  for  the  absorption  at  each  sight  line,  the 
minimum  covering  factor  for  this  line  is  then  25%.  It  is  therefore 
possible  that  the  OH  responsible  for  the  65  /tm  absorption  does 
not  entirely  coincide  with  that  producing  the  84  //m  absorption 
but  is  only  a  fraction  of  the  latter,  consistent  with  its  lower  energy 
level  being  at  nearly  300  K.  Nevertheless,  for  the  sake  of  sim¬ 
plicity,  we  assume  in  this  section  that  both  lines  arise  in  the  same 
region — one  that,  on  the  basis  of  the  84  /im  OH  strength,  covers 
the  total  FIR  continuum.  The  derived  OH  column  densities  will 
be  lower  limits,  and  the  inferred  properties  of  the  continuum 
source  will  be  associated  with  at  least  ~50%  of  the  observed  FIR 
emission. 

The  equivalent  widths  W  are  then  given  by 

ur  1  XI  i/{T Qx)0  I" 

W  =  2  1 - — -  I  (I  -  exp  {-Tv})dv,  (2) 

where  B,,(Tex )  is  the  blackbody  emission  at  the  excitation  temper¬ 
ature  Tex  of  the  line,  O  =  7rR2q/£>2  is  the  solid  angle  subtended 
by  the  source,  F \  is  the  observed  continuum  flux  density  at  the 

2  The  effective  width  is  defined  here  as  /( 1  —  exp  {—  r„})rft>,  where  r„  — 
Toexp{—  (o/AF)2}  and  To  is  the  line  opacity  at  line  center. 
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Fig.  5. — Predicted  equivalent  widths  (W)  vs.  the  OH  excitation  temperature  (Tex)  for  the  119  (solid  curves),  84  (dotted  curves),  65  (dash-dotted  curves),  and  53  (dashed 
curves)  /im  OH  doublets.  Results  are  shown  for  three  assumed  sizes  of  the  continuum  emission,  corresponding  to  the  continuum  models  A,  B,  and  C  in  Table  2.  The  observed 
values  are  indicated  with  horizontal  lines.  An(OH)  is  the  assumed  OH  column  density  in  units  of  1017  cm-2.  The  equivalent  widths  are  positive  if  the  lines  are  in  absorption;  a 
drop  in  W  with  increasing  Tex  thus  indicates  that  the  corresponding  line  is  turning  from  absorption  to  emission. 


wavelength  X  of  the  line,  r„  is  the  line  opacity  at  velocity  v,  and 
the  factor  2  accounts  for  the  two  lines  that  compose  a  doublet. 
The  values  of  W  are  positive  for  absorption  lines,  and  negative 
for  lines  observed  in  emission.  Equation  (2)  applies  both  to  op¬ 
tically  thin  and  optically  thick  lines.  For  optically  thin  lines,  W  is 
proportional  to  the  assumed  column  density;  for  very  optically 
thick  lines,  IFbecomes  insensitive  to  the  column  density  and  scales 
linearly  with  the  turbulent  velocity.  Based  on  DS98,  we  have 
adopted  A  V  =  40  km  s  “ 1 ,  which  is  probably  accurate  within  50%. 

The  fractional  level  populations  and  column  densities  (and  hence 
the  opacities),  and  blackbody  temperatures  required  to  obtain  W 
from  equation  (2)  can  be  estimated  from  a  given  column  density 
A(OH),  and  by  assuming  that  the  excitation  temperature  Tex 
is  the  same  for  all  OH  transitions.  This  assumption  is  certainly 
only  an  approximation.  The  OH  levels  are  pumped  primarily 
through  absorption  of  FIR  photons;  this  is  quite  a  general  model 
result  when  both  the  ground  and  the  excited  OH  rotational  lines 
(except  the  163  pm  one)  are  observed  in  absorption  (see  the  sim¬ 
ilar  conclusion  in  Paper  I  for  the  case  of  Arp  220,  as  well  as  the 
case  of  the  H20  lines  at  25-45  pm  in  Orion/IRc2;  Wright  et  al. 
2000).  When  radiative  excitation  dominates,  the  level  popula¬ 
tions  tend  toward  equilibrium  with  the  radiation  field  within  the 
inner  cloud,  and  the  excitation  temperatures  are  similar  for  all  tran¬ 
sitions  in  those  regions.  However,  the  observed  absorptions  are 
produced  close  to  the  cloud  boundary  where  the  OH  molecules  are 
illuminated  from  only  one  side  and  Tex  is  less  than  the  radiation 
temperature.  Due  to  trapping  effects,  the  lines  with  the  highest 
opacity  (119  and  84  pm  OH  lines)  remain  more  excited  than  the 
other  lines  at  inner  locations  of  the  region,  but  their  Tex  values 
decrease  more  steeply  outwards  and  fall  below  the  Tex  values  of 
other,  thinner  lines  close  to  the  cloud  boundary.  The  detailed,  non¬ 
local  radiative  transfer  models  described  in  §  3.3  show  that  the 
assumption  of  equal  Tex  is  only  approximately  valid  for  the  119 
and  84  pm  lines,  but  for  the  other  lines,  and  in  particular  for  the 
65  and  53  pm  lines,  Tex  may  be  higher  or  lower  depending  on 
position  in  the  region.  Nevertheless,  the  opacities  of  the  119,  84, 
65,  and  53  pm  lines  are  mostly  determined  by  the  Tex  of  the  119 
and  84  pm  lines.  As  a  result,  the  following  analysis  of  the  ra¬ 


diation  density  (or,  equivalently,  Req,  required  to  account  for  the 
observed  absorptions)  is  accurate,  at  least  to  a  first  approximation. 

The  three  parameters  now  required  to  estimate  W  from  equa¬ 
tion  (2)  are  then  A(OH),  Req,  and  Tex.  Figure  5  shows  the  expected 
values  of  W formodels  A,  B,  and  C  as  a  function  of  Tex  forthe  119, 
84,  65,  and  53  pm  lines,  and  compares  them  with  the  observed 
values.  Each  panel  assumes  a  value  of  Req  that  corresponds  to  the 
continuum  models  A,  B,  and  C  described  in  §  3. 1  (see  also  Table  2). 
The  values  of  A'TOH )  used  for  each  model  are  just  reference  val¬ 
ues  discussed  below. 

The  spectral  line  analysis  of  Model  A  (Fig.  5a)  shows  that  this 
scenario  can  be  ruled  out  as  the  main  source  of  FIR  radiation 
from  Mrk  23 1  (Fig.  4a):  the  continuum  model  predicts  low 
(between  20  and  34  K  for  the  different  lines)  but  the  observed 
absorption  in  the  65  pm  line,  with  lower  level  energy  of  290  K, 
requires  a  much  higher  Tex.  Furthermore,  the  119  and  84  ;im  lines 
are  expected  to  be  in  emission  (W  <  0)  as  soon  as  Tex  becomes 
higher  than  20  and  25  K  (which,  on  the  other  hand,  is  not  possible 
for  radiative  excitation).  With  such  a  low  value  of  Tex,  is 

negligible  even  with  our  adopted  screen  A(OH)  =  10 18  cm~2,  a 
value  that  overestimates  (F53/Jm  by  more  than  a  factor  of  2. 

The  above  problems  still  remain  to  some  extent  in  model  B, 
when  Req  is  reduced  to  400  pc  (Fig.  5b).  Here  the  dust  radia¬ 
tion  temperatures  allow  the  lines  to  be  seen  in  absorption  up  to 
Tex  «  40  K;  however,  W65/im  is  underestimated  by  more  than  a 
factor  of  2  for  the  adopted  A(OH)  =  3  x  10 17  cm-2,  yet  this 
column  density  still  overestimates  the  absorption  of  the  53  pm  line. 
Although  model  B  cannot  account  for  the  65  pm  line  strength,  a 
region  of  similar  size  but  lower  A(OH)  could  contribute  to  the 
observed  absorptions  of  the  119,  84,  and  53  pm  lines. 

The  single-component  model  that  best  accounts  for  the  four 
observed  OH  lines  is  model  C,  with  Req  =  200  pc  (Fig.  5c).  The 
corresponding  continuum  model  (Fig.  4c),  with  Tj  =  1A  K,  also 
fits  the  overall  FIR  continuum  emission  rather  well.  Significantly, 
our  models  in  §  3.3  show  that  the  excitation  temperatures  required 
to  reproduce  the  observed  equivalent  widths,  40-60  K,  are  those 
computed  at  the  cloud  surface  if  the  OH  is  excited  by  the  infra¬ 
red  emission  from  a  blackbody  at  Tj  =  74  K.  Finally,  the  dust 
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Fig.  6. — Continuum-subtracted  spectmm  of  Mrk  231  compared  with  two  models  for  the  OH,  H2O,  [C  n],  and  [0 1]  lines.  The  models  are  the  single-component  C  (black 
dashed  line),  and  the  composite  model  D  ( gray  solid  line )  whose  parameters  are  listed  in  Tables  2  and  3 ,  and  whose  predicted  continuum  spectral  energy  distributions  are  shown 
in  Figs.  4c  and  4 d.  The  lines  that  most  contribute  to  the  modeled  spectra  are  labeled. 


temperature  and  gas  mass  ( Table  2)  in  model  C  are  consistent  with 
the  gas  temperature  and  H2  mass  derived  by  P1W07  from  the  sub¬ 
millimeter  CO  and  HCN  emission.  They  found  that  this  warm 
gas  component  hosts  most  of  the  molecular  mass  in  the  Galaxy. 
The  H2  column  density,  IV(H2)  ~  1.5  x  1024  cm~2,  indicates 
high  optical  depths,  as  in  the  Galactic  Sgr  B2  molecular  cloud, 
but  Mrk  23 1  is  much  warmer.  If  the  column  density  in  Mrk  23 1  is 
concentrated  in  a  face-on  disk  of  thickness  H  —  23  pc,  as  con¬ 
cluded  by  DS98,  the  expected  density  is  n(  H2)  ~  2  x  104  cm~3, 
just  the  amount  needed  to  account  for  the  CO  submillimeter  lines 
(P1W07).  On  the  other  hand,  if  this  warm  and  dense  component 
is  identified  with  the  inner  disk  of  radius  460  pc  reported  by  DS98, 
the  area  filling  factor  is/  ~  0.2.  In  spite  of  the  general  agreement 
between  our  model  C  with  other  observations,  a  closer  inspection 
of  this  model  (§  3.3)  reveals  some  discrepancies  with  other  OH 
and  H20  lines  that  suggest  that  a  slightly  modified  scenario  can 
better  explain  the  overall  observed  absorption  patterns. 

3.3.  Models  for  OH  and  H20 

Radiative  transfer  modeling  of  the  observed  OH  and  H20 
lines  was  done  using  the  code  described  in  Gonzalez- Alfonso  & 
Cemicharo  (1997,  1999),  which  computes  the  statistical  equi¬ 
librium  populations  of  a  given  molecule  in  spherical  symmetry. 
Line  broadening  is  assumed  to  be  caused  by  microturbulence. 
Our  code  accounts  for  a  nonlocal  treatment  of  the  radiative  trap¬ 


ping  in  the  molecular  lines  and  of  the  excitation  through  absorp¬ 
tion  of  photons  emitted  by  dust,  as  well  as  for  collisional  excitation. 
Both  line  and  continuum  opacities  for  photons  emitted  in  both 
lines  and  continuum  are  taken  into  account.  Collisional  rates  were 
taken  from  Offer  et  al.  (1994)  and  Green  et  al.  (1993)  for  OH  and 
H20,  respectively.  As  we  also  found  for  Arp  220  (Paper  I),  the 
overall  excitation  is  dominated  by  absorption  of  FIR  continuum 
photons  in  all  models.  If  shock  conditions  (high  density  and  tem¬ 
perature)  were  assumed,  only  the  absorption  in  the  lowest-lying 
lines  would  be  significantly  affected.  Once  the  continuum  model 
is  fixed,  our  results  only  depend  on  the  molecular  column  den¬ 
sities  and  turbulent  velocity  (see  Paper  1  for  a  fuller  description). 

As  mentioned  above  (§  3.2),  the  observed  absoiption  strengths 
are  not  sensitive  to  the  amounts  of  OH  and  H20  in  the  inner 
regions  of  the  modeled  regions,  but  only  to  the  amounts  of  OH 
and  H20  that  are  close  to  the  cloud  (or  disk)  boundary.  For  this 
reason,  we  calculate  two  values  for  the  derived  molecular  column 
densities:  Ascr(X)  denotes  the  column  density  for  a  shell  of  spe¬ 
cies  X  covering  the  infrared  source  (i.e.,  the  screen  case),  whereas 
Nmtx(X)  is  the  inferred  column  density  for  models  where  X 
and  dust  are  evenly  mixed  (the  mixed  case).  Evidently  Amlx  will 
be  much  higher  than  NSCI,  but  from  our  data  there  are  only  a 
few,  nondefinitive  ways  to  discriminate  between  the  alternatives. 
The  163  fin i  OH  and  120  //m  18OH  lines  are  stronger  in  the 
mixed  case,  but  neither  of  these  features  is  unambiguously  detected. 
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TABLE  3 

Models  for  Molecules  and  Atoms 


D  D 

Model3  C  (Warm)  (Cold) 


A(H2)  (cm-3) .  1.5  x  1024  7.2  x  1023  6.5  x  1023 

AV  (km  s-1) .  40  60  40 

Ascr(OH)b  (cm-2) .  1.0  X  1017  9.2  x  1016  6.1  x  1016 

Amix(OH)b  (cm-2) .  4.5  X  1018  1.4  x  1018  9.7  x  1017 

X(OH)c .  3.0  x  10-6  2.0  x  10-6  1.5  x  10-6 

Ascr(H20)  (cm  2) .  5.4  X  1016  6.1  x  1016  <5.0  x  1015 

Nmix( H20)  (cm-2) .  2.0  x  1018  9.6  x  1017  <1.0  x  1017 

X(H20)c .  1.3  x  10-6  1.3  x  10-6  <1.5  X  10-7 

Nmis(C+)d  (cm-2) .  2.3  x  102°  3.1  x  1019e  2.8  x  1019e 

Amix(0)d  (cm-2) .  4.9  x  102°  6.7  x  1019e  6.0  x  1019e 


3  Models  correspond  to  the  continuum  models  listed  in  Table  2  and  shown 
in  Figs.  4c  and  4 d. 

b  (Vscr(X)  is  the  column  density  of  a  shell  of  species  X  surrounding  the  con¬ 
tinuum  source,  and  Nma(X)  is  the  column  density  for  the  case  that  species  X  is 
evenly  mixed  with  the  dust. 

c  Abundances  are  given  relative  to  H2  using  the  .Vnm(  X )  values. 

d  The  column  densities  of  C+  and  0  are  computed  by  assuming  gas  at  400  K, 
and  a  0°  to  C+  abundance  ratio  of  2.1. 

e  Equal  averaged  abundances  are  assumed  in  the  warm  and  cold  compo¬ 
nents  of  model  D. 

Nevertheless,  we  do  not  find  any  strong  arguments  for  thinking 
that  OH  and  H20  are  only  present  on  the  surface  of  the  disk,  and 
so  the  Nmtx  values  may  be  considered  somewhat  more  reliable. 
The  abundances  we  derive  below  are  based  on  this  assumption;  we 
revisit  the  “mixed”  case  when  we  discuss  models  for  the  [C  n]  line. 

Since  model  C  (Fig.  4c,  Fig.  5c,  Table  2)  gives  the  best  single¬ 
component  fit  to  most  of  the  OH  equivalent  widths,  we  first  check 
if  it  can  account  for  the  observed  OH  and  H20  absorption  features. 
Figure  6  compares  the  observed  continuum-subtracted  spectrum 
and  the  modeled  results  (dashed  spectrum,  mixed  case)  for  the 
wavelength  ranges  where  the  signal-to-noise  ratio  is  adequate. 
Table  3  lists  the  physical  parameters  obtained  for  this  model.  We 
have  assumed  a  turbulent  velocity  A  V  of  40kms_1  (§  3.2;DS98). 
The  model  fits  satisfactorily  the  OH  119,  84,  65,  and  53  /tm  lines, 
thus  demonstrating  the  approximate  validity  of  the  simple  method 
outlined  in  §  3.2.  The  value  of  Ascr(OH)  =  1017  cm~2  is  also 
the  same  as  estimated  from  the  equivalent  widths.  The  H20  col¬ 
umn  densities  are  determined  by  the  strengths  of  the  330  — >  22i 
and  331  — >  2  2o  lines  at  66-67  ;im. 

Some  features  of  model  C,  however,  are  inconsistent  with  the 
data.  The  possible  emission  in  the  OH  1 63  /tm  line  is  not  repro¬ 
duced,  and  the  absorption  in  the  79  and  99  /tm  OH  lines  appears 
excessive.  The  model  also  predicts  too  much  absorption  in  the 
322  ->•  2 1 1  (90  /tm),  2 20  ->•  I11  (101  /t: m),  221  -»•  ho  (108  /tm), 
and  4 14  — *  3  03  (113  /tm)  H20  lines.  All  these  discrepancies  sug¬ 
gest  that  the  component  that  accounts  for  the  absorption  of  the 
65-68  /tm  OH  and  H20  lines  is  weaker  than  postulated  in  model 
C  at  wavelengths  longer  than  80  /tm. 

These  discrepancies  may  be  resolved  by  invoking  two  dif¬ 
ferent  components  for  the  FIR  continuum  emission,  as  in  model 
D  (Fig.  4,  Table  2).  The  warm  compact  component,  responsible 
for  the  65-68  /tm  OH  and  H20  lines,  will  produce  weak  absorp¬ 
tions  in  the  80-120  /tm  range  as  a  consequence  of  the  relatively 
weak  continuum  emission  at  these  wavelengths.  The  more  ex¬ 
tended  component  will  contribute  to  the  observed  absorptions  in 
the  53,  84,  and  119  /tm  OH  lines.  The  compactness  of  the  warm 
component  suggests  that  it  is  relatively  close  to  the  AGN,  and 
thus  we  have  assumed  AV  =  60  km  s_1  for  Dwaim  (Table  3); 
this  is  the  turbulent  velocity  found  by  DS98  around  the  rota¬ 


tion  curve  turnover  radius  of  75  pc.  For  the  extended  component 
(Dcoid),  A  V  =  40  km  s  “ 1  is  assumed.  Figure  6  shows  that  a  better 
fit  to  the  overall  spectrum  is  indeed  found  with  this  composite 
model  ( gray  line),  with  the  lines  in  the  80-120  /tm  range  brought 
down  to  levels  compatible  with  observations.  In  addition,  the 
model  predicts  the  163  /tm  OH  line  to  be  in  emission.  The  col¬ 
umn  densities  in  model  D  are  significantly  lower  than  in  C  and,  to 
avoid  the  above  mentioned  H20  absorptions  at  90-120  /tm,  only 
an  upper  limit  for  the  H20  column  density  has  been  derived  for 
the  extended  component.  On  the  basis  of  this  improved  fit  to  the 
molecular  lines,  we  favor  model  D  over  model  C. 

3.4.  Models  for  [C  11]  and  [O  1] 

We  now  check  whether  our  preferred  models  for  OH  and  H20 
are  consistent  with  the  observed  [C  11]  158  /tm  and  [O  1]  63  /tm 
lines,  and  in  particular  if  the  version  with  OH  and  dust  coexis¬ 
tent  (the  mixed  case),  which  predicts  high  OH  column  densities 
( Table  3),  is  compatible  with  the  low  [C  n]  to  FIR  flux  ratio  seen 
and  previously  reported  in  Mrk  231  (Luhman  et  al.  2003).  In 
Paper  1  we  assumed  a  C+  to  OH  abundance  ratio  of  100  (based 
on  standard  gas-phase  carbon  abundance,  and  the  abundance  of 
OH  found  in  the  galactic  center  and  Sgr  B2),  and  found  a  remark¬ 
able  agreement  between  the  predicted  and  observed  [C  11]  line 
flux  in  Aip  220.  In  the  present  case  we  have  searched  for  the  re¬ 
quired  C+  and  O0  column  densities  in  models  C  and  D  (mixed 
case)  that  fit  the  observed  [C  11]  and  [0 1]  line  fluxes,  with  the  ad¬ 
ditional  constraint  that  the  [O  1]  145  /tm  line  is  not  detected  to 
a  limit  of  4  x  10~21  W  cm~2  (see  also  Luhman  et  al.  2003).  In 
model  D,  the  C+  abundance  was  assumed  to  be  the  same  in  both 
components.  Both  the  [C  11]  and  [O  1]  lines  are  expected  to  form 
within  A  y  ~  3  mag  of  the  surfaces  of  PDRs,  and  so  we  allowed 
the  gas  temperature  to  vary  between  250  and  700  K  (e.g.,  Kaufman 
et  al.  1999);  the  results  for  the  [C  11]  line  are  not  sensitive  to  7*  as 
long  as  7)-  is  sufficiently  higher  than  the  energy  of  the  upper  level, 
«90  K.  In  the  regions  directly  exposed  to  the  far-UV  radiation 
field,  all  of  the  carbon  and  oxygen  are  expected  to  be  in  their  ion¬ 
ized  and  atomic  forms,  respectively.  Therefore  an  O0  to  C+  abun¬ 
dance  ratio  of  2.1,  corresponding  to  O0  and  C+  gas-phase 
abundances  of  3  x  10~4  and  1.4  x  1CU4  (Sembach  &  Savage 
1996),  is  imposed.  Using  the  collisional  rates  by  Launay  &  Roueff 
(1977a,  1977b)  for  carbon  and  oxygen,  we  then  searched  for  the 
density  »(H)  and  column  densities  Nmix(C+)  and  /Vmlx(O0)  re¬ 
quired  to  match  the  line  strengths.  The  results  for  both  models  C 
and  D  are  also  in  Figure  6,  and  the  calculated  column  densities 
are  listed  in  Table  3 . 

We  could  not  find  any  satisfactory  fit  to  the  three  atomic  lines 
with  model  C.  Since  the  continuum  is  optically  thick  at  X  <  150  /tm, 
extinction  effects  are  important  in  the  [0 1]  63  /tm  line,  and  densities 
higher  than  5  x  103  cm~3  were  required  to  account  for  it.  On  the 
other  hand,  high  column  densities  were  also  needed  to  fit  the 
[C  11]  line,  and  the  final  result  is  that  all  combinations  that  fitted 
both  the  [C 11]  and  the  [0 1]  63  /tm  lines  yielded  more  flux  than  the 
upper  limit  for  the  undetected  [O  1]  145  /mi  line.  The  model  in 
Figure  6  shows  the  result  for  h(H)  =  8  x  103  cnr3  and  7*  = 
400  K.  Only  high  densities  (~105  cnr3)  and  relatively  low  tem¬ 
peratures  (Si  150  K)  can  approximately  match  both  O0  lines,  but 
then  the  [C  11]  line  is  underestimated  by  more  than  a  factor  of  2. 

In  contrast,  model  D  allowed  us  to  find  a  more  satisfactory  fit 
to  the  three  lines.  As  we  argue  below,  the  compact-warm  compo¬ 
nent  is  expected  to  give  negligible  emission  in  the  C+  and  O0 
lines.  The  extended  component  Dco\a,  with  its  moderate  dust  tem¬ 
perature  ( Td  =  47  K)  and  continuum  opacity,  can  explain  the 
[0 1]  63  /un  line  with  a  density  and  column  density  lower  than  in 
model  C,  thereby  predicting  a  weak  145  /tm  line.  The  model  in 


312 


GONZALEZ-ALFONSO  ET  AL. 


Vol.  675 


Figure  6  uses  »(H)  =  5  x  1 0 3  cnr3  and  7*  =  400  K  to r  /7>co|tl.  It 
is  remarkable  that  when  we  generated  the  same  model  but  excluded 
the  effect  of  dust,  we  found  that  the  [O  i]  63  /tm  line  was  stronger 
by  a  factor  of  4.2.  The  dust  has  a  small  effect  on  the  other  two 
[O  i]  145  /tm  and  [C  n]  158  /tm  lines,  but  can  substantially  lower 
the  [O  i]  63  /mi  line  because  the  absorption  of  63  /mi  line-emitted 
photons  by  dust  grains  (i.e.,  extinction)  is  more  important  than 
the  additional  pumping  of  0°  atoms  by  63  /tm  dust-emitted 
photons.  Therefore,  the  ratio  of  the  [C  n]  158  /tm  line  to  the 
[O  i]  63  /im  line,  observed  of  order  unity,  is  reproduced  with 
higher  densities  and  temperatures  than  in  previous  studies: 
Figure  4  of  Kaufman  et  al.  (1999)  shows  that  7([0  i]  63  /mi)/ 
/([C  n]  158  /tm)  «  1  for  »(H)  =  5  x  103  cnr3  and  Go  ~  1 02, 
corresponding  to  a  temperature  somewhat  higher  than  100  K 
(Fig.  1  of  Kaufman  et  al.  1999).  Including  dust  leads  to  a  sim¬ 
ilar  intensity  ratio  for  7*-  -  400  K,  or  Go  <  10 '.  On  the  other  hand, 
our  model  for77coid  yields  /([ C  n]  158  /tm)//([0  i]  145  /tm)  ss 
7.5,  which  implies  Go  ~  103  for  the  above  density  (see  Fig.  6  of 
Luhman  et  al.  2003).  These  results  suggest  that  the  effect  of  ex¬ 
tinction  cannot  be  neglected  when  modeling  the  [O  i]  63  /tm  line 
in  sources  that  are  optically  thick  at  wavelengths  shorter  than 
100  /tm.  Self-absorption  by  foreground  low-excitation  clouds 
may  also  lower  the  flux  of  the  [O  i]  63  /tm  line  and  to  a  lesser 
extent  that  of  the  [C  n]  line  (Fischer  et  al.  1 999;  Vastel  et  al.  2002; 
Luhman  et  al.  2003;  Paper  I). 

To  compare  our  results  for  line  intensity  ratios  involving  the 
[O  i]  63  /tm  line  with  the  corresponding  plots  shown  by  Luhman 
et  al.  (2003)  one  should  use  the  versions  of  our  models  that  ig¬ 
nore  dust,  because  the  intensities  given  in  the  models  by  Kaufman 
et  al.  (1999)  are  uncorrected  for  extinction.  In  our  models  that 
ignore  dust,  [C  n]/[0  i]  ~  0.4  and([C  n]  +  [O  i])/Ifir  ~  2  x  10~3. 
Here,  the  FIR  flux  corresponds  to  that  of  component  Dcoid 
(Table  2).  Using  these  values  in  Figure  4  of  Luhman  et  al.  (2003) 
one  obtains  «(H)  =  a  few  x  103  cm~3  and  Go  =  a  few  x  103, 
roughly  consistent  with  our  input  model  parameters. 

The  C+  and  O0  column  densities  required  to  fit  the  lines  in 
model  77 cold  are  rather  high.  The  averaged  C+  and  O0  abundances 
are  2.2  x  10~5  and  4.6  x  10~5  relative  to  H  nuclei,  respectively. 
If  the  C+  to  OH  abundance  ratio  is  100  (Paper  I),  the  mixed  col¬ 
umn  density  of  OH  is  overestimated  by  a  factor  of  3 .  In  any  case, 
the  Nmm{C+)  values  indicate  that  an  important  fraction  of  C 
nuclei  is  in  ionized  form,  ^  15%  in  model  77coid  for  a  gas  phase 
carbon  abundance  of  1 .4  x  10~4  (Sembach  &  Savage  1996).  The 
corresponding  mass  of  atomic  gas  is  8  x  108  Me,  at  least  1  order 
of  magnitude  higher  than  the  amount  estimated  in  M82  and  NGC 
278  (Wolfire  et  al.  1990;  Kaufman  et  al.  1999).  As  mentioned 
above,  C+  is  expected  to  be  abundant  only  in  regions  within  A 
2  mag  from  the  surfaces  of  PDRs,  and  the  relatively  high  frac¬ 
tion  of  C+  derived  above  indicates  that  those  regions  must  be 
widespread  in  the  Mrk  23 1  disk.  On  the  other  hand,  one  possible 
interpretation  of  the  high  H2  column  densities  required  to  gen¬ 
erate  the  radiation  field  needed  to  account  for  the  OH  lines  is  that 
we  are  viewing  a  number  of  PDRs  that  overlap  along  the  line  of 
sight  (see  also  §  4). 

We  now  check  the  possibility  that  a  significant  fraction  of  the 
[C  n]  158  /tm  emission  arises  from  a  more  extended  region,  such 
as  the  outer  disk  found  by  DS98  in  CO.  For  an  outer  disk  mass  of 
2  x  1O9M0  and  density  of  102  cm~3  (DS98),  and  using  an  upper 
limit  for  the  C+  abundance  of  <10  4  and  collisional  rates  at 
250  K  (from  Launay  &  Roueff  1977a),  we  calculate  that  5/1/4 
of  the  observed  [C  ii]  line  arises  from  this  extended  region.  There¬ 
fore  we  do  not  expect  much  contribution  from  it,  unless  the  C+  is 
concentrated  in  clumps  of  much  higher  density. 


Our  model  results  can  be  used  to  re-examine  the  issue  of  the 
[C  n]  line  deficit  found  in  ULIRGs.  Contrary  to  what  might  be 
expected  from  the  low  [C  n]  to  FIR  luminosity  ratio  in  Mrk  23 1 , 
relatively  high  average  abundances  of  C+  are  required  in  /Jc,m  to 
reproduce  the  [C  n]  line.  Explanations  based  on  self-absorption 
effects,  extinction  by  dust,  effects  of  dust-bounded  H  n  regions, 
high  Go/m,  and  high  formation  rate  of  CO  at  the  expense  of  C+ 
have  been  invoked  to  account  for  the  [C  n]  deficit  (Luhman  et  al. 
1998, 2003;  Malhotra  et  al.  2001;  PIW07;  Paper  I,  and  references 
therein).  In  our  models  for  Mrk  231,  the  high  abundance  of  C+ 
yields  the  observed  low  [C  n]  to  FIR  luminosity  ratio  because  of 
the  high  luminosity-to-gas  mass  ratio  (L/M)  of  our  continuum 
models.  In  both  model  C  and  model  D,  an  overall  L/M  ss  500 
Lq/Mq  is  derived  (Table  2),  where  the  gas  mass  includes  both 
the  atomic  and  molecular  components.  We  argue  that,  if  this  ex¬ 
treme  L/M  ratio  applies  to  the  whole  region  from  which  the  FIR 
emission  arises,  as  appears  to  be  the  case  for  Mrk  23 1 ,  the  [C  n] 
line  strength  relative  to  the  infrared  emission  can  never  attain 
values  in  excess  of  0.1%.  The  [C  n]  luminosity  emitted  per  unit 
of  total  gas  mass  is  given  by 


</>[Cn]  =  ,  ,  ,  7  v  — nr  ,nc+Aulhv,  (3) 
1  +  (9i/9u)  exP  {Eu/kTex} 

where  cj  l  and  gu  are  the  degeneracies  of  the  lower  and  upper  level, 
Eu  is  the  upper  level  energy,  Au/  is  the  Einstein  coefficient  for 
spontaneous  emission,  and  nc+  denotes  the  number  of  C+  ions 
per  unit  of  total  gas  mass  (including  all  components).  It  is  as¬ 
sumed  in  equation  (3)  that  the  line  is  optically  thin.  In  the  limit 
T’ex  »  Eu/k  one  obtains 


0[Cu]  ~  °-5 


X(C+)  Lp 
10-4  A70  ’ 


(4) 


where  y(C+)  denotes  the  average  C+  abundance  relative  to  total 
H  nuclei.  Therefore,  the  expected  [C  n]-to-IR  emission  is 


^[Cn]  _  3x(C  +  ) 
L/M  ~  10-4 


(5) 


for  L/M  =  500  Le/Me. 3  Hence,  even  assuming  the  highest  nearly 
possible  x(C+)  ~  10~4,  the  [C  n]-to-lR  ratio  is  not  expected  to 
be  above  10~3  for  such  high  values  of  L/M.  The  reason  is  simply 
the  limited  reservoir  of  gas,  and  hence  of  C+  ions,  per  unit  of  out¬ 
put  power.  High  L/M  values  are  in  Mrk  23 1  primarily  due  to  the 
contribution  of  the  AGN,  but  could  also  be  associated  with  high 
star  formation  efficiencies  (maybe  extreme;  see  Scoville  2004): 
if  we  attribute  the  total  infrared  luminosity  to  FUV  heated  grains 
in  PDRs  and  use  the  units  of  Go  (i.e.,  the  Habing  Field  of  1 .6  x 
1 0  3  ergs  cm~2  s  1 )  to  calculate  the  L/M  for  a  PDR  column  den¬ 
sity  of  Av=  10  mag  with  the  standard  gas-to-extinction  conver¬ 
sion  A ( H  +  2H2)/H„  ~  2.5  x  1021  cm~2  mag-1  (e.g.,  Dickman 
1978),  we  find  that  the  value  of  Go  that  corresponds  to  L/M  = 
500  Lq/Mq  is  2.5  x  1 04,  which  we  express  for  the  general  case 
as 


—  =  500 
M 


Go 


2.5  xlO4 


10  mag 


W’  <61 


3  For  simplicity,  the  luminosity  accounts  here  for  the  total  infrared  emission 
as  derived  in  our  models,  whereas  in  previous  studies  the  FIR  emission  was  de¬ 
fined  between  40  and  120  /im;  the  correction  factor  is  less  than  2  (see  Sanders  & 
Mirabel  1996). 
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so  that  high  L/M  values  would  indicate  high  Go  values.  For  PDRs 
with  Go  ~  2.5  x  104  and  Av  ~  10  mag,  the  average  C+  abun¬ 
dance  is  ~2  xlO-5  (i.e.,  carbon  is  ionized  within  ~2  mag  of  the 
surfaces  of  PDRs)  and  <j>  [c„i  /(L/M)  will  attain  the  lowest  values 
observed  in  ULIRGs. 

In  our  favored  model  D,  the  component  Dwarm  (L/M  =  3270 
Lq/Mq),  which  emits  most  of  the  IR  luminosity,  could  only  ac¬ 
count  for  the  [C  n]  line  if  essentially  all  gas-phase  carbon  were  in 
ionized  form.  However,  this  seems  implausible,  as  the  C+  abun¬ 
dance  in  XDRs  is  y(C  '  )  ~  a  few  x  1 0  5  (Meijerink  &  Spaans 
2005),  and  FUV  photons  from  stars  are  absorbed  in  thin  shells 
of  dust  on  the  surfaces  of  PDRs.  With  the  assumed  column  den¬ 
sities  in  Table  3,  the  calculated  contribution  from  Dwarm  to  the 
observed  [C  n]  emission  is  only  8%,  the  [O  i]  145  /imflux  is  15% 
of  the  upper  limit,  and  the  [O  i]  63  pm  is  predicted  in  weak  ab¬ 
sorption  against  the  continuum.  ( In  this  calculation,  7*  =  103  K 
and  n  =  2  x  103  cm~3  are  assumed,  so  that  the  atomic  region 
is  in  pressure  equilibrium  with  the  cooler  molecular  gas  with 
Tk  ~  102  Kand«  =  2  x  1 04  cm  3.)  It  is  therefore  expected  that 
the  strongest  source  of  infrared  radiation  (the  warm  component) 
gives  low  [C  n]  and  [O  i]  line  emission  even  if  one  assumes 
relatively  high  averaged  fractional  abundances  for  C+  and  O0.  As 
in  Luhman  et  al.  (2003),  we  then  invoke  a  component  that  emits 
strongly  in  the  FIR  but  yields  weak  emission  in  the  [C  n]  and  [O  i] 
lines.  In  contrast  to  Luhman  et  al.  (2003)  however,  we  leave  open 
the  possibility  that  Dwarm  has  a  relatively  high  C+  abundance,  as 
expected  from  both  XDR  and  PDR  chemistry.  In  starburst  and 
normal  galaxies,  the  moderate  star  formation  efficiencies  (e.g., 
~15Z,0/M0  inNGC  1068;  Planesas  etal.  1991),  implying  lower 
values  of  Go,  will  not  place  such  a  strong  limit  on  the  C+  cooling 
rate  relative  to  the  infrared  emission.  We  suggest,  then,  that  the 
[C  n]  158  pm  line  deficit  found  in  ULIRGs  may  be  associated 
with  a  high  L/M  ratio  that,  allowing  a  relatively  high  C+  abun¬ 
dance,  indicates  an  important  contribution  to  the  FIR  emission 
by  an  AGN,  and/or  a  high  star  formation  efficiency. 

4.  DISCUSSION  AND  CONCLUSIONS 

Simultaneous  modeling  of  the  FIR  continuum  and  line  absorp¬ 
tions  is  the  only  consistent  way  to  derive  from  FIR  observations 
the  physical  and  chemical  properties  of  bright  infrared  galax¬ 
ies,  where  due  to  high  radiation  densities,  the  OH  and  H20  lines 
are  pumped  through  absorption  of  continuum  photons.  In  Mrk  231, 
the  high-lying  65-68  pm  OH  and  H20  lines  indicate  the  pres¬ 
ence  of  a  compact  (^100  pc)  and  warm  (^100  K)  region,  pre¬ 
sumably  located  around  the  AGN,  with  high  continuum  opacity 
(t ioo/im  ~  1).  The  size  of  this  region  is  comparable  to  (although 
somewhat  larger  than)  the  inner  nuclear  disk  where  the  OH  mega¬ 
masers  have  been  observed  from  VLBI  observations  (Klockner 
et  al.  2003),  and  is  similar  to  the  spatial  extent  of  the  regions  where 
the  bulk  of  the  OH  megamaser  emission  is  generated  (regions  C, 
NE,  and  SW  in  Fig.  5  of  Richards  et  al.  2005).  This  strongly 
suggests  that  we  are  observing  the  FIR  counterpart  of  the  18  cm 
OH  emission  (e.g.,  Skinner  et  al.  1997).  Richards  et  al.  (2005) 
estimate  A(OH)  ~  1.5  x  10 15  cm~2  in  1  pc  clouds  where  AV  « 
1.7  km  s_1;  therefore,  A(OH)/ A V  ~  1015  cnr2  (km  s-1)-1.  In 
both  the  warm  and  cold  components  of  model  D,  we  have  ob¬ 
tained  NSCT(OH)/AV  ~  1.5  x  1015  cm~2  (km  s-1)-1,  in  rough 
agreement  with  the  estimates  by  Richards  et  al.  (2005).  The 
“screen”  case  is  most  appropriate  for  comparison  here  because 
the  OH  megamasers  amplify  the  background  radio  continuum 
emission. 

The  nuclear  region  is  surrounded  by  a  more  extended  (Req  ~ 
350  pc)  and  colder  (~45  K)  region,  which  produces  the  bulk  of 
the  [C  n]  and  [O  i]  line  emissions  and  contributes  to  the  low-lying 


OH  lines.  Allowing  for  an  area  filling  factor  of  0.6-0. 7,  the  to¬ 
tal  extent  of  the  FIR  emission  is  then  comparable  to  the  size  of 
the  CO  (inner)  disk  ( Biyant  &  Scoville  1 996;  DS98).  The  inferred 
densities  [m(H2)  ~  104  cm~3]  may  explain  the  submillimeter 
CO  emission  detected  by  PIW07  and  somewhat  lower  densities 
(a  few  x  103  cm~3)  can  account  for  the  [O  i]  line  emission.  We 
do  not  however  rule  out  the  possibility  that  some  of  the  FIR  emis¬ 
sion  arises  from  a  more  extended  region,  such  as  the  outermost 
disk  observed  by  DS98  in  the  millimeter  CO  lines  or  the  much 
more  extended  region  observed  in  soft  X-rays  (Gallagher  et  al. 
2002). 

The  FIR  emission  from  Dwarm  is  probably  dominated  by  the 
AGN,  as  pointed  out  in  §  3.1.  Davies  et  al.  (2004)  estimated 
that,  within  a  radius  of  330  pc,  the  starburst  luminosity  accounts 
for  (0.8- 1.3)  x  10 12  L0;  since  the  luminosity  for  the  component 
Dcoid  is  «7.5  x  1011  Le,  the  starburst  luminosity  contribution 
fromDwaral  is  most  probably  ;£6  x  10 11  Le.  Davies  et  al.  (2004) 
also  estimated  that,  within  100  pc  of  the  AGN  (the  size  ofDwarm), 
the  starburst  contribution  is  (3. 2-4. 7)  x  1011  L0,  in  agreement 
with  the  afore  mentioned  upper  limit.  On  the  other  hand,  the  ex¬ 
tended  (360x260  pc2)  radio  continuum  1.4  GHz  emission  de¬ 
tected  by  Carilli  et  al.  (1998)  is  compatible  with  both  the  whole 
FIR  emission  arising  from  the  starburst  (yielding  an  IR-to-radio 
flux  density  ratio  of  Q  ~  2.5;  Carilli  etal.  1998),  and  only  the  Dcoid 
FIR  emission  arising  the  starburst  (yielding  Q  =  2.1).  If  one 
now  assumes  that  the  starburst  contribution  to  the  FIR  emission 
ofDwarm  is  3  x  10 11  Le,  the  resulting  Q  is  2.3  (in  agreement  with 
the  median  value  of  Q  =  2.3  ±  0.2  in  normal  galaxies;  Condon 
1992,  and  references  therein).  In  addition,  the  surface  brightness 
and  luminosity-to-mass  ratio  of  Dwarm  due  to  only  the  starburst 
component  would  then  be  2  x  10 12  Le  kpc~2  and  550  Lq/Mq, 
which  are  probably  upper  limits  for  a  starburst.  In  summary,  our 
preferred  model  D  favors  that  ~2/3  of  the  bolometric  luminosity 
from  Mrk  23 1  is  due  to  the  AGN. 

Although  the  chemistry  of  the  extended  starburst  Dcoid  com¬ 
ponent,  where  the  lower-lying  OH  lines  and  the  bulk  of  the  [C  n] 
158  pm  and  [O  i]  63  pm  lines  are  generated,  appears  to  be  dom¬ 
inated  by  PDRs,  the  dominant  chemistry  of  the  warm  component 
may  be  of  a  different  nature.  Since  the  compactness  of  Dwarm  sug¬ 
gests  that  its  FIR  emission  is  dominated  by  the  AGN,  and  given 
the  high  X-ray  intrinsic  luminosity  of  the  AGN  (~1044  erg  s_1) 
in  the  2-10  keV  band  (Braito  et  al.  2004),  we  consider  the  pos¬ 
sibility  that  the  derived  OH  and  H20  column  densities  can  be 
better  explained  in  terms  of  XDR  chemistry.  Calculations  by 
Meijerink  &  Spaans  (2005)  show  that,  in  dense  XDRs  with  Ah  ~ 
1024  cm~2,  the  OH  and  H20  column  densities  are  expected  to 
attain  values  of  ^1018  cm~2,  similar  to  the  values  that  we  have 
derived  in  the  mixed  case.  Furthermore,  the  Compton-thick  screen 
with  Ah  ~  2  x  1024  cm~2  required  to  block  the  primary  X-ray 
emission  (Braito  et  al.  2004)  could  be  partially  identified  with 
Dwarm,  with  similar  column  density.  Therefore  an  XDR  may  nat¬ 
urally  explain  the  observations  of  high-lying  OH  and  H20  lines 
in  Mrk  23 1 .  Nevertheless,  both  OH  and  H20  are  also  tracers  of 
starburst  chemistry,  and  a  number  of  PDRs  and  hot  cores  over¬ 
lapping  along  the  line  of  sight  could  in  principle  yield  the  col¬ 
umn  densities  required  to  explain  our  observations  at  least  in  the 
screen  case,  and  could  then  mimic  an  XDR.  We  conclude  that 
further  observations  are  still  required  to  discern  between  XDR  and 
starburst  chemistry. 

We  suggest  that  the  [C  n]  158  pm  line  deficit  in  Mrk  23 1 ,  and 
also  probably  in  Arp  220  (see  Paper  I),  is  primarily  due  to  a 
component  that  dominates  the  FIR  emission  but  emits  weakly 
in  the  [C  n]  line,  as  previously  suggested  by  Luhman  et  al.  (2003). 
However,  we  also  suggest  that  this  component  may  still  be  rich  in 
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C+,  but  with  an  extreme  luminosity-to-gas  mass  ratio  that  limits 
the  [C  n]  luminosity  per  unit  of  luminous  power  in  the  continuum. 
Our  derived  luminosity  to  H2  mass  ratio  of  500  Le/M@  is  high; 
future  studies  will  indicate  whether  it  can  be  applied  to  other 
ULIRGs  with  similar  [C  n]  deficits. 

Given  the  high  C+  column  densities  (a  few  x  10 19  cm~2)  re¬ 
quired  to  account  for  the  [C  n]  158  /mi  line  emission,  the  H2 
columns  derived  toward  the  extended  component  are  unlikely 
to  be  associated  with  non-overlapping  star- forming  regions  sur¬ 
rounded  by  optically  thick  envelopes.  A  single  PDR  is  not  ex¬ 
pected  to  have  N(C+)  in  excess  of  a  few  x  10 17  cm~2.  A  scenario 
with  a  crowded  population  of  PDRs  overlapping  along  the  line 
of  sight  is  therefore  more  plausible.  If  a  typical  single  OB  stel¬ 
lar  cluster  emits  up  to  10 39  erg  s  the  maximum  value  found  by 
Scoville  et  al.  (2001)  in  M51  and  also  the  luminosity  of  W49, 
then  the  number  of  single  clusters  to  produce  the  Mrk  23 1  star- 
burst  luminosity  of  10 12  Le  is  ^4  x  106.  If  these  clusters  are 
concentrated  in  a  disk  of  radius  460  pc  and  thickness  25  pc  (DS98), 
the  mean  distance  between  neighbor  clusters  is  only  ^  1 .5  pc  (see 
also  Keto  et  al.  1992).  Merging  of  PDRs  may  be  compatible  with 
some  area  filling  factor  if,  for  example,  some  spiral  structure  is 
invoked  within  the  disk. 

In  NGC  1068,  Spinoglio  et  al.  (2005)  observed  the  OH  119  //m 
line  in  emission  against  the  continuum.  On  the  basis  of  a  possible 
XDR  chemistry,  and  given  the  corresponding  mass,  density,  and 
temperature  that  characterize  the  nuclear  region,  Spinoglio  et  al. 
(2005)  suggested  that  the  OH  119  /tm  emission  line  could  be 
formed  in  that  region.  This  option  also  relied  on  the  fact  that  rel¬ 
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atively  weak  FIR  emission  is  expected  to  arise  from  the  circum- 
nuclear  disk  of  the  Seyfert  2  galaxy,  so  that  the  OH  119  /tm  line 
can  be  excited  through  collisions  and  emit  above  the  continuum 
in  these  dense  and  warm  environments.  The  situation  is  differ¬ 
ent  for  Mrk  23 1  and  Arp  220.  In  these  ULIRGs,  the  extremely 
high  luminosity  arising  from  the  nuclear  region,  together  with 
the  high  concentrations  of  gas  there,  make  the  nuclear  FIR  emis¬ 
sion  component  very  bright.  Any  possible  emission  in  the  1 1 9  /tm 
line  will  be  obscured  by  the  strong  absorption,  which  in  these 
objects  is  ^30  times  stronger  than  the  emission  feature  in  NGC 
1 068.  In  addition,  the  strong  FIR  radiation  density  pumps  higher- 
lying  OH  and  H20  levels,  thus  producing  absorptions  in  the 
higher-excitation  lines.  The  OH  and  H20  molecules  are  there¬ 
fore  potentially  powerful  tracers  of  circumnuclear  regions  around 
AGNs.  Future  Herschel  observations  of  Mrk  23 1  and  other  sources 
will  allow  us  to  apply  the  models  developed  in  this  paper  to  these 
sources,  refine  them  accordingly,  and  will  certainly  give  new  in¬ 
sights  into  the  physical  and  chemical  conditions  of  bright  in¬ 
frared  galaxies. 
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